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Abstract 

We study the formation of high energy tails in a one-dimensional kinetic model 
for granular gases, the so-called Inelastic Maxwell Model. We introduce a time- 
discretized version of the stochastic process, and show that continuous time implies 
larger fluctuations of the particles energies. This is due to a statistical relation 
between the number of inelastic collisions undergone by a particle and its average 
energy. This feature is responsible for the high energy tails in the model, as shown 
by computer simulations and by analytical calculations on a linear Lorentz model. 
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1 Introduction 



Lorentz gas models have been used historically in order to clarify the fea- 
tures of more complex non-linear kinetic equations. For instance, P. and T. 
Ehrenfest (jl|) introduced the so-called wind tree model, in order to discuss 
the Stosszahlansatz of the Boltzmann equation. In this paper, we will use a 
similar approach in order to study inelastic gases, namely low-density systems 
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composed by a large number of moving macroscopic particles, themselves per- 
forming inelastic interactions. Due to their inelasticity, inelastic gases dissipate 
kinetic energy. This implies that the granular temperature, defined kinetically 
by T~<y^>, where V = v — u(r, t) and u(r, t) is the local mean velocity, 
asymptotically vanishes if the system is not supplied by an external energy 
source. Before reaching total rest state, these systems usually reach a self- 
similar solution, i.e. form preserving solution whose time dependence occurs 
through the granular temperature, /(f ; t) = —^=fs{—^=). Such scaling ve- 
locity distributions have been observed in a large variety of kinetic models, 
and have been shown to genericallv highlight non-Maxwellian features and 
overpopulated high energy tail (0; 0)- Let us note that the specific shape of 
the tail may exhibit very different behaviours, depending on the details of the 
model (j^. 

In this paper, we show on simplified kinetic models that such high energy tails 
may originate from the fact that the average energy of particles depends on 
their collision history, i.e. on the number of collisions they have undergone in 
the course of time. This effect, that has no counterpart in the case of elastic 
interactions, is verified by Direct Simulation Monte-Carlo (DSMC) simulations 
of the non-linear Inelastic Maxwell Models, and by focusing on a discrete 
time version of the dynamics. Analytical calculations are also performed for a 
simpler linear Lorentz model. 



2 Inelastic Maxwell Model 



In the following, we study scaling solutions in the context of the one-dimensional 
Inelastic Maxwell Model (IMM). This model (0) derives from a mathematical 
simplification of the non-linear Boltzmann equation for inelastic hard rods (6:), 
assuming that the collision rate may be replaced by an average rate, indepen- 
dent of the relative velocity of the colliding particles. We also assume that the 
system is and remain homogeneous in the course of time, so that the resulting 
kinetic equation writes: 



^^(^i; ^) + /(^^; t)=l dv, {-)f{v[- t)f{v,- 1) 



dt J a 



where the primed velocities are the pre-coUisional velocities, defined by the 
collision rule: 

(1 + a) / (1 + a) 

Vi = Vi Vi2 f2 = ^2 H 7, ^12 (2) 

Za la 



2 



Vi2 is the relative velocity vi — V2. The restitution parameter a E ]0 : 1] mea- 
sures dissipation of energy at each collision, and the elastic limit corresponds 
to the case a = 1. It is well-known that the scaling solution of 1 reads (|7|): 



/T(t)' 



where the granular temperature decreases like T{t) = Tq e 



Recently, we have introduced (|80 a time discretized version of 1: 

oo ^ 

fN+livi) = I dv2 i-)fNiv[)fNiv'^) (4) 

J a 

— oo 



Let us call continuous time dynamics (CTD) and discrete time dynamics 
(DTD) the models associated to 1 and 4 respectively. Arguments based on 
the central limit theorem show that any initial velocity distribution converges 
towards a scaling solution whose shape is gaussian in the case of DTD: 

fNiv) = (5) 

(27rTAr) 2 



Tjv is the temperature after N collision steps, Tjy = (^-^)^To. One should 
stress that the transition from a power law tail 3 towards a Maxwell-Boltzmann 
distribution 5 is entirely due to the passage from a continuous time towards 
a discrete time description. The main difference between these two models 
follows. In DTD, every particle has suffered exactly inelastic collisions at 
time A^, so that it has dissipated the same fraction of energy on the average 
at that time. In contrast, in CDT, the constituents of the inelastic gas have 
performed different number of inelastic collisions at a given time t, dis- 

J.N —t 

tributed according to the Poisson law . This suggests that ensembles of 
particles, discriminated by their number of collisions, should be characterised 
by different quantities of energy on average, thereby increasing fluctuations of 
the particles energies. 

In order to verify this possible explanation, we have performed DSMC simu- 
lations (jH) of the kinetic equations 1 and 4, based on their interpretation as 
a stochastic process. The CDT case is well-known and consists in picking a 
given number of random collision pairs at each step. The DTD case is based on 
the same algorithm, except that at each step, the N particles are distributed 
randomly into y collision couples. This method ensures that each particle 
collides one and only one time at each step. At the end of each step, the veloc- 
ities are rescaled in order to keep the total energy constant Simulation 
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Fig. 1. In (a), velocity distributions in arithmetic and in logarithmic scale (inset) for 
DTD. The solid line corresponds to the theoretical Maxwell-Boltzmann distribution. 
In (b), the simulations of CDT show the convergence towards the theoretical solution 
3 (solid line). In both cases, the system is composed by 150000 particles with a=0.5. 
The asymptotic velocity distribution is plotted with circles (o). The dashed lines 
are guides for the eye in order to isolate the initial conditions, that are uniform or 
Maxwell-Boltzmann distribution. 

methods have been verified by looking at the asymptotic velocity distribution, 
and checking the theoretical results 3 and 5 (fig.l). In the case of CDT, we 
have also highlighted a non-trivial relation between the number of collisions 
suffered by a particle, and its average energy. In order to evaluate this relation, 
we have started a simulation from a Maxwell-Boltzmann initial condition, and 
let the simulation run during 50 collisions per particle. After that time, the 
simulation is stopped, the particles velocities are rescaled so that the average 
energy < E >= 1 and we measure the average energy < E >k of particles 
having performed K collisions. The results (fig. 2) clearly show that the aver- 
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Fig. 2. Average energy of particles as a function of their number of collisions. The 
dashed line is the constant value < E >= 1. The system is composed of 20000000 
particles with a = 0.5. 

age energy of particles is a decreasing function of their number of collisions, 
thereby confirming the above discussion. 



3 One dimensional Lorentz model 



Lorentz-hke models have already been applied to the study of inelastic gases, in 
order to show non-equipartition of energy for instance (|l8|). These systems are 
composed by particles that do not interact between themselves, and undergo 
collisions with randomly distributed static (infinitely heavy) scatterers. In this 
paper, we will consider a model which is very similar to the Ehrenfest wind 
tree model (jl|): 

^^ = \v\i-,f{v';t)-fiv;t)) (6) 
ot 



where the pre-coUisional is defined by v = —^v. This kinetic equation has 
been shown (2; 3) to be related to the inelastic Liouville equation (14) and to 
population dynamics In the following, we use an IMM-like approximation 
in order to treat this kinetic equation, i.e. we assume that the collision rate 
|f I may be replaced by a mean value 



'T(t) and we rescale the time scale. We 
call the resulting model the Inelastic Lorentz Walk (ILW): 



df{v;t) 
dt 



-f{-;t)-f{v;t) 
a a 



d(f){k;t) 

at 



[ak; t) — (f){k; t) 



(7) 
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The right side equation is the kinetic equation for the characteristic function 
of the velocity distribution (f){k;t) = ^ J dvf{y,t)e^^'". In Fourier space, the 
discrete time version of the process reads: 

(t>N+i{k) = (t>N{oik) (8) 
whose solution is: 

Mk) = Mc^^'k) (9) 

Let us stress that the stochastic process leaves the solution scale invariant, i.e. 
9 is a self-similar solution whose form is that of the initial velocity distribu- 
tion and whose temperature asymptotically vanishes. Obviously, the discrete 
process does not lead to the formation of high energy tails, but this property 
ceases to be true if only a fraction p of the discs performs a collision at each 
step. In this case, the state of the system after one step writes: 

Mk) = il-P) Mk)+pM(^k) (10) 

Simple analytical calculations show that the kurtosis k, = ^^2^2 of the velocity 
distribution increases through 10 whatever the initial velocity distribution and 
proportion 1. This is due to the fact that the superposition of identical 
distributions with different temperatures Ti, T2 has a higher energy tail than 
the same distribution with the mean temperature r=(l-p) Ti + p To ( lot 17). 



Contrary to the non-linear IMM, the link between the discrete and the con- 
tinuous time dynamics is straightforward. Indeed, usual methods of random 
walk theory (jlj) lead to: 

00 fN -t 

mt) = T.^uo^''k) (11) 



The system at time t is thus composed by different classes of particles that 
are characterised by their mean energy a^, or equivalently by the number 
of undergone collisions. The proportion of particles in these classes is given by 
the probability for a particle of having performed collisions before time t, 
namely ^, . On should note that 11 is a superposition of the initial distri- 
bution at different energies, as in 10 and leads to the formation of high energy 
tails by the same mechanism. Assuming that the small k development of the 
initial characteristic function is (po{k;t) = I]^o formal solution 11 

reads: 

mt) = T.-^e-('--'^^ (12) 
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Detailed analysis of this solution shows that the inelastic Lorentz walk exhibits 
multiscaling properties, i.e. the dimensionalized velocity moments rrii = (< 
{v'^y >)~i decay with different cooling rates nii ~ e~^^^. The cooling rates 
are given by Xi = ^—f- and are decreasing function of i: Aj < Xj, i > j. 
Consequently, the system is characterised by an infinite number of independent 
cooling rates. Moreover, the velocity moments rrii, i > 2, grow towards infinity 
as compared to mi. 



4 Conclusions 

In this paper, we study the formation of the high energy tails observed in 
inelastic gases. To do so, we focus on the one-dimensional Inelastic Maxwell 
Model, which is a mean field approximation of the Boltzmann equation for 
inelastic hard rods. By comparing a discrete and a continuous time version of 
the process, we show how fiuctuations of the number of collisions imply the 
emergence of high energy tails. To do so, we perform DSMC simulations of 
both dynamics, thereby highlighting a new relation between number of inelas- 
tic collisions undergone by a particle and its average energy. This relation 
is also analytically studied by focusing on a simpler linear model, i.e. the In- 
elastic Lorentz Walk. Let us stress that this mechanism is specific to inelastic 
gases, where energy is dissipated at each collision, and has no counterpart in 
elastic gases. A generalisation of this work to higher dimensional systems and 
to more general kinetic models is under progress. 
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